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Abstract. We prove that the electromagnetic material parameters are 
^ I uniquely determined by boundary measurements for the time-harmonic 

. Maxwell equations in certain anisotropic settings. We give a uniqueness 

' result in the inverse problem for Maxwell equations on an admissible 

, Riemannian manifold, and a uniqueness result for Maxwell equations 

^ ' in Euclidean space with admissible matrix coefficients. The proofs are 

based on a new Fourier analytic construction of complex geometrical 
optics solutions on admissible manifolds, and involve a proper notion of 
' uniqueness for such solutions. 

(N ■ 



1. Introduction 



<, 

(— I I Let {M,g) be a compact Riemannian manifold with smooth boundary 

dM, and assume that dimM = 3. We consider the inverse problem of 
recovering electromagnetic material parameters of the medium (M, g) by 
probing with time-harmonic electromagnetic fields. The fields in (M, g) are 
described by complex 1-forms E and H (electric and magnetic fields), and 
, the behavior of the fields is governed by the Maxwell equations in M, 

^ : ( *dE =iu;fiH, 

^ ; \ *dH = -iujeE. ^ 

i Here u; > is a fixed frequency, d is the exterior derivative, and * is the 

' Hodge star operator on {M,g). The material parameters are given by the 

complex functions e and (permittivity and permeability, respectively) . We 
I assume the following conditions on the parameters: 

e,^GC~(M), (1.2) 

^ '. Re(e) > 0, Re(/x) > in M. (1.3) 

■ For the inverse problem, we need to describe the electromagnetic field 

measurements at the boundary dM. Let i : dM — > M be the canonical 
embedding, and consider the tangential trace on A;- forms, 

t : n''{M) n''{dM), 1] ^ i*r]. 

There is a discrete set of resonant frequencies such that if to is outside this 
set, then for any / in Q^{dM) the system (II. ID has a unique solution {E, H) 
satisfying tE = f (see Theorem lA.ip . We shall assume that 

a; > is not a resonant frequency. (1-4) 
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The boundary measurements are given by the admittance map 



The inverse problem for time-harmonic Maxwell equations is to recover the 
material parameters e and n from the knowledge of the admittance map A. 

In the case of lossy materials, one writes e = Re(e) + ia/uj where o" > 
is the conductivity. The zero frequency case (that is, w = 0) then formally 
corresponds to the conductivity equation 



Here 6 is the codifferential. In three and higher dimensions, the inverse 
problem of determining a from boundary measurements for the conductivity 
equation was studied in [4] in a special class of Riemannian manifolds. 

Definition. A compact 3-manifold {M,g) with smooth boundary dM is 
called admissible if (M, g) is embedded in (T, g) where T = R x Mq, (Mq, (70) 
is a simple 2-manifold, and g = c{e(Bgo) where c is a smooth positive function 
and e is the Euclidean metric on R. 

Simple manifolds are defined as follows: 

Definition. A compact manifold (Mo,go) with smooth boundary OMq is 
called simple if for each p in Mq the map exp^ is a diffeomorphism from a 
closed neighborhood of in TpM^ onto Mq, and if OAIq is strictly convex 
(meaning that the second fundamental form of OMq is positive definite). 

Admissible manifolds include compact submanifolds of Euclidean space, 
hyperbolic space, and minus a point, and also sufficiently small subman- 
ifolds of conformally flat manifolds. If M is a bounded open set in R'^ with 
smooth boundary, equipped with a metric which in some local coordinates 
X = {xi,x') has the form 



then {M,g) is admissible if gQ is a simple metric in some sufficiently large 
ball. Also, admissible manifolds are stable under small perturbations of gQ. 
See [1] for more details. 

We will prove the following result, showing that boundary measurements 
for the Maxwell equations uniquely determine the material parameters in 
an admissible manifold. 

Theorem 1.1. Let {M,g) be an admissible manifold, and let (ei,/ii) and 
{^2,fJ'2) be two sets of coefficients satisfying (ll.2|) ~ (ll.4|) . If the admittance 
maps satisfy Ai = A2, then £1 = £2 and /ii = ^2 in 

The second result involves Maxwell equations in a bounded domain Vt in 
R^ with smooth boundary. The coefficients e, // are assumed to be smooth 
positive definite symmetric (1, l)-tensors. Associated to these tensors are 
traveltime metrics g^ and g'^, which are Riemannian metrics in describing 
propagation of waves with different polarizations. We shall assume that the 
velocity of wave propagation is independent of polarization, which amounts 
to the property that e and // are in the same conformal class |12j . 



A : Q:^{dM) n 



(dM), tE^tH. 



6{adu) = 0. 
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The Maxweh equations in 17 can be written as 

( 'V X E = iuJuH, ,^ . 

\ V X ^ = -iujeE, ^ ■ ' 

where E and H are complex vector fields and > is a fixed frequency. 
We consider the electric boundary condition 

-E'tanbn = /, (1-6) 

where / is a smooth tangential vector field on and -Etanlan is the tan- 
gential part of E\qq^. Under the above assumptions, there is a discrete set of 
resonant frequencies outside which the boundary problem for Maxwell equa- 
tions has a unique smooth solution {E,H) (see Section [7]). The admittance 
map is given by 

A : Eti,n\dn ^ Hta,n\dn- 
The next result considers the inverse problem of recovering the electromag- 
netic parameters from A. 

Theorem 1.2. Let £j and fxj be smooth symmetric positive definite (1, 1)- 
tensors on 0,, and suppose that > is not a resonant frequency for the 
corresponding boundary problems. Let Aj be the corresponding admittance 
maps {j = 1,2). Assume that there is a fixed admissible metric g m il. such 
that El, /xi, and H2 are conformal multiples of g^^. If the admittance 
maps satisfy Ai = A2, then £1 = £2 and ^1 = ^2 in ^■ 

To our knowledge, Theorems 11.11 and 11.21 are the first positive results 
on the inverse problem for time-harmonic Maxwell equations in anisotropic 
settings. For bounded domains in where g is the Euclidean metric. 
Theorems 11.11 and 11.21 were proved in |17j . 

There has recently been considerable interest in invisibility cloaking [7], 
where one looks for anisotropic materials for which uniqueness does not hold. 
The prescriptions of electromagnetic parameters for cloaking [6j satisfy that 
e = fi. Moreover the parameters are singular, so that one of the eigenvalues 
is zero at the boundary of the cloaked region. Theorems 11.11 and 11.21 imply 
that there is no cloaking for materials whose electromagnetic parameters 
satisfy the given conditions. 

Formally, the proofs of Theorems 11.11 and 11.21 follow the Euclidean case. 
The proof of the uniqueness result in [l7j was considerably simplified in [19] , 
and the simplified proof can be described by the following seven steps: 

1. Reduction of the Maxwell system to a Dirac system, by introducing 
two auxiliary scalar fields $ and ^. A solution X of the Dirac system 
gives a solution to the original Maxwell system iff $ = ^ = 0. 

2. Reduction to a rescaled Dirac system {P — k + W)Y = 0, where Y 
is obtained by rescaling the components of X by e^^^ and fi^^'^. 

3. Reduction to the Schrodinger equation (—A — k"^ + Q)Z = 0, which 
is possible since {P - k + W){P + k - W^) = -A-k"^ + Q. 

4. Construction of complex geometrical optics solutions to the equation 
(- A -k'^ + Q)Z = 0, which also gives solutions Y = {P + k- W^)Z 
to the Dirac system. 
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5. Construction of solutions to the original Maxwell system. This re- 
quires showing that the scalar fields in Step 1 vanish identically, 
which follows from a uniqueness result for Z. 

6. Inserting complex geometrical optics solutions in an integral identity, 
which allows to recover nonlinear differential expressions involving 
the electromagnetic parameters. 

7. An application of the unique continuation principle for a semilinear 
elliptic system to recover the parameters. 

In |18j . it was shown that Steps 1 to 3 above can be carried out also for 
the Maxwell equations on a Riemannian manifold {M,g). However, Step 4 
requires complex geometrical optics solutions, and these were only available 
for the Euclidean metric. Therefore, it was not possible to go further in the 
non-Euclidean case. 

A construction of complex geometrical optics solutions for scalar elliptic 
equations, valid on admissible Riemannian manifolds (M,g), was given in 
[1]. We will combine the ideas in [1] with the scheme outlined above to 
prove the uniqueness result for the inverse problem for Maxwell equations 
on admissible manifolds. 

It will turn out that the main technical obstacle is Step 5, which requires 
a uniqueness result for the complex geometrical optics solutions. In [1] the 
construction of solutions is based on Carleman estimates, and there is no 
concept of uniqueness for the solutions so obtained. In this article we give 
a new construction of solutions based on direct Fourier arguments. This 
construction comes with a suitable uniqueness result, which can be used to 
carry out the proof of the Maxwell result. 

The main step in the new construction is a counterpart of the basic norm 
estimates of Sylvester-Uhlmann |22j . We outline the idea in a simple case. 
The estimate is valid in {T,g) where T = R x Mq and g = e (B go, but here 
(Mo,5o) can be any compact (n — l)-dimensional manifold with boundary 
(no restrictions on the metric). We look for solutions of the equation 



with Ag the Laplace-Beltrami operator in (T, g) and r a large parameter. 

In the Sylvester-Uhlmann estimates T = R" and g is the Euclidean 
metric, / is in a weighted space such that {x)^^^f G L^(R"') where 
— l<d<0, and one obtains a unique solution u with {x)^u € L^(RP). 
Here 



In our case we write xi for the special Euclidean coordinate in T, and use 
Agmon-type weighted spaces 



and H}q{T) is the set {u £ H}{T) ; uIrxSMq = 0}. 

The next result is a special case of Proposition 14.11 (since there is no 
potential it follows that one may take tq = 1). 




(1.7) 
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Theorem 1.3. Let 6 > 1/2. If |t| > 1 is outside a discrete set, then for any 
/ G L'g{T) there is a unique solution u £ H\^q{T) of the equation (jl.7p . In 
fact, one has u G H'^g{T) and 

MH^_^(T)<C\Tr'\\f\\qiT), 0<S<2, 

with C independent of r. 

In the Sylvester-Uhlmann result, the proof applies the Fourier transform 
and one obtains uniqueness by fixing decay at infinity. In our case there is a 
transversal metric in Mq, and the Fourier transform or conditions at infinity 
are not readily available. However, one can ask for decay in the Euclidean 
variable and Dirichlet boundary values on OMq. This makes it possible to 
use the Fourier transform in xi and eigenfunction expansions in Mq. 

The proof of Theorem 11.31 is robust in the sense that one can essentially 
replace the Laplacian in Mq by any positive operator with a complete set 
of eigenfunctions. We will need this flexibility in the Maxwell result when 
proving similar estimates for the Hodge Laplacian on forms. There is also 
an extra twist in the construction of solutions since one needs a result like 
Theorem 11.31 which applies to functions / which may not decay (so one is 
out of the standard Agmon setting) , see Sections H] and [5] for these more 
general results. 

The construction could also be used to develop constructive methods for 
certain anisotropic inverse problems. In the Euclidean case, results of this 
type were given in [16] for the 3D conductivity equation and in [T7] for 
Maxwell equations. 

Earlier work on the inverse problem for the Maxwell system in Euclidean 
space includes a study of the linearized inverse problem [20j , a local unique- 
ness result [21], and a result for the corresponding inverse scattering prob- 
lem in the case where fi is constant [2]. As mentioned above, the full inverse 
problem was solved in [TT] , and in [TO] the proof was simplified and also a 
reconstruction from measurements based on dipole point sources was given. 
The paper [IH] is a survey and also considers the manifold setting. The in- 
verse problem for Maxwell equations in chiral media was considered in [T^] . 
Boundary determination results are given in |10J and [14]. Finally, [l] gives 
a partial data result for this problem, based on Isakov's method [9]. For 
results on inverse problems for the Maxwell equations in time domain, we 
refer to [12] and the references therein. 

The structure of the paper is as follows. Section [2] contains notation 
and identities in Riemannian geometry which will be used throughout the 
article. The reductions of the Maxwell equations to Dirac and Schrodinger 
equations are given in Section [3l The norm estimates and uniqueness results 
required for constructing complex geometrical optics solutions are given in 
Sections HI and [S] and the construction of solutions is taken up in Section [Gj 
In Section [7| we prove Theorems 11.11 and 11.21 There are two appendices, one 
on the wellposedness theory of boundary value problems for Maxwell, and 
one including a unique continuation result for principally diagonal systems 
required for our results. 
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2. Notation and identities 

We will briefly introduce some basic notation and identities in Riemannian 
geometry which will be used throughout. We refer to [23] for these facts. 

In this section let (M, g) be a smooth (= C°°) n-dimensional Riemannian 
manifold with or without boundary. All manifolds will be assumed to be 
oriented. We write {v,w) for the (7-inner product of tangent vectors, and 
\v\ = for the g-norm. If x = (xi, . . . local coordinates and 

dj the corresponding vector fields, we write gjk = {dj,dk) for the metric in 
these coordinates. The determinant of {gjk) is denoted by \g\, and {g^^) is 
the matrix inverse of (gjk)- 

We shall often do computations in normal coordinates. These are coor- 
dinates X defined in a neighborhood of a point p £ M™* such that x{p) = 
and geodesies through p correspond to rays through the origin in the x 
coordinates. The metric in these coordinates satisfies 

fl'jfc(O) = 5jk, digjk{0) = 0. 

For points p S dM we will employ boundary normal coordinates, which are 
coordinates y = {y',yn) near p so that y{p) = 0, y' are normal coordinates 
on dM centered at p, and ynig) is the geodesic distance from a point q to 
dM. The metric has the form 

5(y) = (^°^^^ 1), go{y) = {9My))Tk=v 

and gjk{0) = 5jk, digjk{0) = 0. We denote by u the 1-form corresponding to 
the outer unit normal vector of dM, so that = —dy"^ in boundary normal 
coordinates. 

The Einstein convention of summing over repeated upper and lower in- 
dices will be used. We convert vector fields to 1-forms and vice versa by the 
musical isomorphisms, which are given by 

(A^a,)^ = Afc(^x^ Xk = gjkX^, 
{ujk dx^ f = io^ dj , J = g^^Uk . 

The set of smooth /c-forms on M is denoted by Q^M, and the graded algebra 
of differential forms is written as 

QM = el^o^'^M. 

The set of fe-forms with or H'^ coefficients are denoted by L'^{Qf'M) and 
H^{n^M), respectively. Here for s S R are the usual Sobolev spaces on 
M. The inner product and norm are extended to forms and more generally 
tensors on M in the usual way. 

Let d : OJ^M fi^+^M be the exterior derivative, and let * : Vt^M 
Qn-kj^ be the Hodge star operator. We introduce the sesquilinear inner 
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product on i7*^M, 

iv\0= I {v,OdV= I r?A*C. 

JM JM 

Here dV = *1 = \g\^/'^ dx^ ■ ■ ■ c?.t" is the volume form. The codifferential 
5 : fi^M — > Q^~^M is defined as the formal adjoint of d in the inner product 
on real valued forms, so that 

{dvlO = ivl^O^ foi' V £ ^''~^M,( e QI^M compactly supported and real. 
These operators satisfy the following relations on A;-forms in M: 

** = [-ifin-k)^ g ^ ^_^^kin-k)-n+k-l 

If ^ is a 1-form then the interior product is the formal adjoint of in the 
inner product on real valued forms, and on A;-forms it has the expression 

The Hodge Laplacian on fc-forms is defined by 

-A = (d + 5f = dd + 5d. 

It satisfies A* = *A. 

The Levi-Civita connection, defined on tensors in M, is denoted by V. 
We will slightly abuse notation and reserve the expression V/ (where / is 
any function) for the metric gradient of /, defined by 

V/ = = g^^djfdk. 
The and norms may be expressed invariantly as 

ll/lbi(M) = ll/llL2(M) + l|V/||i2(M), 
ll/lk2(M) = ll/lki(M) + l|V2/llL2(M)- 

Here of course ||T'||^2(jv^) = (/Afl^^P dV)^^'^ for a tensor T. 

For n = 3, the surface divergence of / G H^{p?-{dM)) is given by 

J^w{f) = {daMf.dS) 

where dS is the volume form on dM. A computation in boundary normal 
coordinates shows that Div(/) = —{v,*du)\QM where u G H^'^^^'^{^^M) is 
any 1-form with tu = f (here ,s > 0). 

Finally, in the case n = 3, we collect a number of identities which will be 
useful for computations. Below let / be a smooth function, a = aj dx^ and 
/3 = f3j dx^ and 7 = 7^ dx^ three 1-forms, rj a /c-form, and an Z-form. For 
the Hodge star one has 

* * rj = ij, 
*(«A*/?) = {a,P), 
*(a A *[/? A 7]) = {a, 7)/? - {a, 
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and the operators d and 6 satisfy 



Si] = {-!)'' *d*r], 



d{7] A C) 

6{a A /?) = 



5a 



-\9\-'/%{\9\'/'9^'a,), 
= dr|^Q + {-lfr]^dQ, 
f5rj +{-!)'' *df A *r], 
i6a)P-{6(3)a - [q^/3«] 



3. Reduction to Schrodinger equation 



In this section we present the reductions of the time-harmonic MaxweU 
system to Dirac and Schrodinger equations, which corresponds to Steps 1 
to 3 in the introduction. This mostly follows [18J and [19j although with 
different notations. We will also give a reduction to the case where the 
coefficients are constant near the boundary, namely, 

e = Eq and fi = fiQ near dM for some constants eo, /^o > 0. (3-1) 

It is well known that the Maxwell system is not elliptic as it is 

written. We perform an elliptization by adding the constituent equations, 
obtained from (jl.ip by applying d* to both equations: 



Adding two equations requires adding two unknowns, which will be the 
scalar fields ^ and ^. The choice for how to couple $ and ^ into the larger 
system obtained from (jl.ip . (|3.2p was motivated in [19] by dimensionality 
arguments. The end result is the following system: 



Here we have written D = -d and a = log e, P = log fi. We will also write 
D* = — i(5 for the formal adjoint of D in the sesquilinear inner product on 
forms. 

We wish to express (j3.3p as an equation for the graded differential form 
X = ^ + E + *H + written in vector notation as 



Note that we have grouped the even and odd degree forms together. This 
will result in a block structure for the equation. Now, taking Hodge star of 
the first and last equations in ()3.3p results in the system 




(3.2) 



D * E + Da A*E -ujfi*^ 
*D'$ + DE -Ljfi* H + *Da A ^ 

D * H + DP A *H - Loe * 
*D^ - DH + *DI3 A<^-uje*E 



0, 
0, 
0, 
0. 



(3.3) 



X = *H\*^ E) . 



{P + V)X = 



(3.4) 
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where P and V are given in matrix notation by 



/ 


D* \ 






*Da A * 




D* D 


, v = 


— UJfl 


*Da A * 


D 




Df3A 


—uje 


\ D D* 


J 




\ DI3A 


—uje 



p 



This is the first Dirac equation we will use. Note that P is just the self- 
adjoint Dirac type operator D + D* on JIM, and that {E, H) solves the 
original Maxwell system (jl.ip ifi' X solves (|3.4p with ^ = = 0. 

For the reduction to a Schrodinger equation, it will be convenient to 
rescale X by 



X 



(3.5) 



where y = ( yo \Y^ )*, and is the /c-form part of F G nM. 
Assuming for the moment, a direct computation using the identities 
in Section [2] shows that (|3.4|) is equivalent with the rescaled Dirac equation 
for Y: 



{P - k + W)Y = 0. 
Here W is the potential, with compact support in M"'-^, given by 



W 



(3.6) 





/ 




*Da A* \ 


1 






*Da A * —DaA 


2 












*DP A * 


J 



where k = uj{£ij,)^^'^, k = w(eoAio)^''^- 

We will also need the potential W^, which is the formal transpose of W 
in the inner product on real valued forms, given by 



-{K-k) + 





/ 


*D(3 A* \ 


1 




*DP A * -Df3A 


2 


DaA 






\ DaA *Da A * 


J 



The adjoint is W* = W^. The following result contains the Schrodinger 
equations, involving the Hodge Laplacian —A = dS + Sd on QM, in a form 
which will be convenient below. 

Lemma 3.1. We have 

{P -k + W){P + k-W^) = -A-k"^ + Q, 
{P + k - W^){P - k + W) = -A - k'^ + Q', 
{P -k + W*){P + k-W) = -A-k'^ + Q, 
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where Q, Q' , and Q are smooth potentials with compact support in M™*, 



Q = k^- 
Q' = k^- 





/ Aa + 


^{da, da) 





• 




2 1 





• 


• 


• 


K + - 

2 














• 


• 


Ap + ^{dp,dp) 









• 


• 





• / 






h{dP,dp) 










2 1 

K 

2 




• 


• 


• 


• 










Aa — ^{da, da) 







l 


• 


• 


• 


• J 



and • denote smooth coefficients. 

Proof. We give the proof of the first identity, the other ones being analogous. 
One has 

{P - k + W){P + k - W^) = -A - k'^ + W{P + k)-{P- k)W^ - WW\ 
The point is to show that the first order term vanishes. We write W as 

W = -{K-k) + ]^W^, 

where Wq acts on a graded form X = Xj^ + X^ , with X^ = X^ + X^ and 
X_ = X^ by 

WqX = {-Da A +iDa)X- + {D(3 A -iD^)X+. 

We will use the identities in Section [2j If is a 0-form then 

{WoP - PWo*)n = ^{Wodu - {d - 5){uDa)) 

= —{—da A du + *da A *du — du Ada + u{5da) — *du A *da) 
= {Aa)u. 
If n is a 1-form we have 

{WqP - PWli)u = ^{Wo{du - 6u) - {d - 5){-Dfi Au + {D(3,u))) 
= -{d/3 Adu + *d/3 A *du - {6u)df3 - df3 A du - d{df3, u) - 6{d/3 An)). 

The identity 6{dP A u) = (— A/?)n — {6u)dp — [V/?, n"]'' and a computation 
in normal coordinates implies that 

{WqP - PWl^)u = 2{V^p)u - {Ap)u. 

Here {V'^(3)u denotes the 1-form corresponding to the vector field {V'^P)uK 
The computation for 2-forms and 3-forms can be reduced to the previous 
cases by noting that if u is a A;-form, then 

{d-5)*u= (-1)^ *{d + 5)u, {t] a -in) *u= (-1)''"^ * (r? A +in)u. 

Thus, if *u is a 2-form then 

{WqP - PWq*) * u 

= *[da Adu + {5u)da — *da A *du — da Adu + d{da, u) + 5{da A u)] 
= *[2(V^a)n - {Aa)u]. 
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Similarly, if *u is a 3-form then 

{WqP - PWo*) *u = (A/3) * u. 
We have P{fu) — fPu = {Df A +iDf)u for a function /, so 

WP - PW' = ^{WoP - PWo) + A +iD^. 

This shows that {P - k + W){P + k - VF*) = -A - k"^ + Q where Q is an 
operator of order 0. Since 

k{W + W^) - WW^ = k^ - + ^k{Wo + W^) - jWqW^ 

where WqW^o-'^ = - [{Da A -iDa)'^ X+ + {DP '^-iDf})'^ and since one has 
{C A —i^)'^u = —{^,,^)u for any A:-form u, we obtain the required expression 
for Q. □ 

The preceding arguments show how to reduce the original Maxwell system 
to Dirac and Schrodinger equations. In the next lemma, which is similar to 
[19\ p. 1135], we give a reduction on the level of boundary measurements: 
if the admittance maps for two Maxwell systems coincide, then one has an 
integral identity involving the potentials Qj and solutions of the Schrodinger 
and Dirac systems. Note that Zi has to be related to a solution for Maxwell, 
but I2 need not be. This flexibility in the choice of Y2 will simplify the 
recovery of coefficients. 

Lemma 3.2. Let (ei,/ii) and {82,^^2) be two sets of coefficients satisfying 
()1.2p - ()1.4|] . and assume that Ai = A2. After replacing {M,g) by a larger 
manifold (which is admissible if {M,g) is), one may assume that 

£1 = £2 = £0 and Hi = iJ-2 = fJ-o near dM for constants eO)/"o > 0, (3.7) 

and one has the identity 

{{Qi-Q2)Zi\Y2) = (3.8) 
for any smooth graded forms Zj,Yj satisfying the following properties: 
{P-k + Wi)Yi = 0, Yi = {P + k-Wl)Zi, 
{P-k + W^)Y2 = 0, 
yo = y3 ^ 0. 



The reduction to the case where ()3.7p holds is a consequence of the next 
boundary determination result. 

Theorem 3.3. Let (M, g) be a compact 3-manifold with smooth boundary, 
and let e and satisfy ()1.2p -( fT~i]) . Given a point on 5M, the admittance 
map A uniquely determines the Taylor series of e and /x at that point in 
boundary normal coordinates. 



Proof. This result was proved in [T^, [Mj in the case where M is a smooth 
domain in and g is the Euclidean metric (the result was for complex e and 
real //, but the same proof works also for complex //). The argument proceeds 
by showing that the admittance map is a pseudodifferential operator on dM, 
and by computing the symbol of A in boundary normal coordinates at a 
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fixed point p G dM. One then proves, by looking at the difference of two 
admittance maps, that the Taylor series of e and n are uniquely determined 
at p. 

Fortunately, if (M, g) is a Riemannian manifold with boundary, the form 
of the metric in boundary normal coordinates is exactly the same as in the 
Euclidean case. This means that the arguments of [TU], [2], which were 
given in boundary normal coordinates of a Euclidean domain, carry over 
without changes to establish Theorem 13.31 for any Riemannian manifold 
(M,g). □ 

Proof of Lemma \3.SX We first establish (13. 8p in the original manifold M. 
Note that if ()3.7p is not satisfied, we may formally take A; = in the preceding 
arguments and then all conclusions remain valid except that Wj and Qj may 
not be compactly supported in M™*. 

Let Zj and Yj be as described, and let Xi be the solution to (j3.4p . with 
potential Vi, corresponding to Yi as in (j3.5p . Since Ai = A2, we can find a 
solution X2 of {P + V2)X2 = with tH2 = tHi and tE2 = tEi on dM. Here 
we write 

Xi = ( *Hi\Q Ei)\ = ( I ^2 )* • 

If Y2 is the solution to (P — A; + W2)Y2 = corresponding to X2 as in (|3.5p . 
then ^(11 — Y2) = since ei = £2 and ^1 = ;U2 on dM by Theorem 13.31 
We wish to argue that 

A (Fi - I2) = 0, - 1^2) = on 9M. (3.9) 

The first part is immediate since u A r] = on dM iff fr/ = 0. For the 
second part we use the surface divergence. The fact that Xi and X2 solve 
the Maxwell equations, together with Theorem 13.31 implies that 

{iy,Hi - H2)\aM = :^(^'/^r^ * '^^1 - * dE2)\dM 

^ -Div(t(^i - E2)) = 0. 



A similar result is true for E1 — E2. This proves (13. 9p since we have *i^A*r/ 
= (i^; ??) for any 1-form rj. 
Let us next prove that 

{{W^-W2)Yi\Y2)=Q. 

We have 

{{Wi - W2)Yi\Y2) = {WiYi\Y2) - {Yi\W^Y2) 
= {Yi\{P - k)Y2) - {{P - k)Y^\Y2) 

= {Yi\{P - k)Y2) - {{P - k){Yi - Y2)\Y2) - {{P - k)Y2\Y2) 
= (FiKP - k)Y2) - (Yi - Y2\{P - k)Y2) - {{P - k)Y2\Y2). 
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In the last step, the boundary term arising from the integration by parts 
vanishes because of (j3.9p . We obtain 

((T^i - W2)Yi\Y2) = {Y2\{P - k)Y2) - {{P - k)Y2\Y2) 

= -iY2\W^Y2) + iW2Y2\Y2) = 0. 

Now ()3.8p will follow if we can prove that 

{{Wi - W2)Yi\Y2) = {{Qi - Q2)Zi\Y2). 
To show this, we recall that Qj in Lemma l3. II has the form 
Qj. = Wj{P + k)-{P- k)Wj - WjWj. 

Then 

{{Wi - W2)Yi\Y2) = {{Wi -W2)iP + k- W\)Zx\Y2) 

= ((Qi + {P- ^)W\)Z^\Y2) - {{Q2 + {P- k)Wi + W2{Wi - Wi))Zi\Y2) 

= ((Qi - ^2)^111^2) + {WiZi\{P - k)Y2) - {W^ZiliP - k)Y2) 

-{{Wi-Wi)Zi\W^Y2). 

Here, we used that Wl = W2 on dM by Theorem 13.31 so there are no 
boundary terms. Now (j3.8|) follows by using the identity —W2Y2 = {P—k)Y2 
in the last term. 

Finally, we show how it is possible to arrange that (j3.7|) holds. The 
definition of admissible manifolds allows to find (upon enlarging (Afo,5o) if 
necessary) a connected admissible manifold {M,g) such that 

M CC M CC T. 

If M is not required to be admissible then any choice M DD M will do. 
By the condition Ai = A2 and by Theorem 13.31 we may extend Sj and 
fij smoothly to M so that ei = 82 and /ii = ^2 in M \ M, Sj and 
have positive real parts in M, and further for some constants eo;/t^o one has 
ei = £2 = £0 and Hi = ^2 = /^o near dM. 

Let now , Yi , I2 be smooth graded forms in M satisfying the conditions 
in the statement of the lemma in M. Since the restrictions to M satisfy the 
same conditions in M, we have (j3.8p in the set M. However, Qi = Q2 in 
M \ M, so (j3.8p remains valid in M. This proves the lemma upon replacing 
M with M. □ 

4. Norm estimates and uniqueness 

In this section let (Mq, go) be a compact {n — l)-dimensional Riemannian 
manifold with smooth boundary, without any restrictions on the metric. 
Consider the cylinder T = R x Mq with metric g = c(e © go), where e is the 
Euclidean metric on R and c is any smooth positive function in T satisfying 

c{xi,x') = 1 when \xi\ is large. 

Here and below, we write xi for the Euclidean coordinate and x' for coor- 
dinates on Mq. The Laplace-Beltrami operators on {T,g) and (Mo,go) are 
denoted by A = and A^' = A^q, respectively. We will use the space 
L^{T) = L^{T,dVg) and the Sobolev spaces H%T). 
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If 5 G R, define the weighted norms 



U t2 



(T) = \\{^i)^'<^\\lHT), 
MhUT) = \\{xi)^u\\Hs^Ty 



Let Li^iT) and Hf{T) be the corresponding spaces. We also consider the 
spaces i?/oc(T) = {n G L\^^{T) ; u G H^{[-R,R] x Mo) for ah i? > 0} and 

Fio(r) = {n G /7i(r) ; n|RxOAfo = 0}, 

^ioc,o(r) = e hUt) ; ^^iRxeAfo = 0} 

The construction of complex geometrical optics solutions in [3] and [TT] 
is based on limiting Carleman weights. It is shown in ^ that the function 
(p{x) = xi is a natural limiting Carleman weight in {T,g). We consider the 
conjugated Helmholtz operator 

The following result gives a norm estimate, corresponding to the Carleman 
estimate in [H Theorem 4.1], and a uniqueness result for this operator. 

Proposition 4.1. Let k > he fixed, let 6 > 1/2, and let g be a potential 
satisfying {xi)'^^q G L°°(T). There exists tq > 1 such that whenever 

> ■'"0 and t"^ + k'^ ^ Spec(— A^,/), 

then for any / G L'j[T) there is a unique solution u G H\^q{T) of the 
equation 

e^^i(-A- A;2 + g)e-^^in = / in T. (4.1) 
Further, u G //^^(T), and the solution satisfies the estimates 

\\u\\HU(T)<C\Tr^\\f\\L2(T), 0<S<2, 

with C independent of r and /. 

For the proof, we first claim that it is enough to consider the case where 
c = 1. To see this, note that if g = eg where g = e(B go, one has the identity 

n + 2 9 N / n — 2 

c 4 (-Ag -k^ + q){c 4 

= (-A^ -k^+ \k'^{l -c)+cq- c'^Ag(c-'^)j )v. (4.2) 

n — 2 

Consequently, u solves ()4.1|) iff t; = c ■* u solves 

e™i (- A^ - fc2 + g)e-™i^; = c^^f. 

Here g = k'^{l — c) + cq — Ag{c ~) is another potential such that 
since c = 1 for |xi| large. This reduction shows that 
Proposition 14.11 will follow from the special case where c = 1. 

Thus, we assume that c = 1 and initially also q = 0. Then g has the form 

1 

go{x') 

and the equation ()4.ip may be written as 

{-dl + 2Tdi-T'^ -k^ - A^,)u = f. (4.3) 
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We are looking for a solution u with u|r,x9Mo = 0- This motivates the 
partial eigenfunction expansions along the transversal manifold: 

oo oo 

u{xi,x') = ^n(xi, 00/(2;'), f{xi,x') = ^f{xi,l)(l)i{x') 
1=0 1=0 

where (pi are the eigenfunctions of — A^,/ on Mq, satisfying —Ax'(j)i = Xi(j)i in 
Mo and (pi\dMo = 0- 

Inserting the expansions of u and / into (j4.3p results in the equations 

i-df + 2Tdi - - + Xi)u{xul) = f{xi,l). (4.4) 

These are second order ODE for the partial Fourier coefficients. To solve 
them, we will use the following simple result on solutions of linear ODE 
involving Agmon type weights. 

Lemma 4.2. If fi = a + ib where a, b are real, a 7^ 0, consider the equation 

u' - fj,u = f in R. (4.5) 

There is a unique solution u = S^f G o5^'(R) for any / S ^'(R). One has 
: L'j{'R} L'j(R) if 5 G R, and also the norm estimates 

\\S^.f\\LjiK) < SII/IIlKr), I«I > 1 and 5 e R, 
\\S,fhl^^n) < ^^ll/llLi(R), a/ and 5 > 1/2. 



The constant C only depends on 6. 

Proof. We take Fourier transforms in ()4.5|) and observe that for / G ^'(R), 
there is a unique solution u = S^f G o5^'(R) given by 



where m(^) = {i^ — fJ-) ^- The condition a 7^ implies that m is a smooth 
function which satisfies 

||mW||Loo < fc!|or(^+^), = 0,1,2,... . 

Thus, for any 5 G R we have mf G if / G , which implies 5^/ G L^. 
If \a\ > 1 then \\mf\\HS < Q|a|-i ||/||^. and \\S^f\\L2 < Cslar^fW^^ 

o o 

We now assume / G L| for (5 > 1/2. If a > 0, the solution to (14. 5p is 
given by 

POO 

S,f[x) = - f{t)e->^('--Ut. 

J X 

This has the estimate 

foo / foo \ 1/2 

\s^f[x)\ < / \f{t)\dt< / {ty^'dt 



<C5\\f\\ 



Li 



since 5 > 1/2. Thus one has ||5u/||^2 < C5||/||£^2 again since 6 > 1/2. A 

— S s 

similar argument gives the result if a < 0. □ 
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Proof of Proposition \4T1\ As argued above, we may assume c = 1. Let us 
also first take q = 0. Then we are looking for solutions to the equation (j4.3p . 

Let < Ai < A2 < • • • be the Dirichlet eigenvalues on —A^' in Mq, and 
let (pi G Hq(Mq) be the corresponding eigenfunctions normalized so that 
{4>i}iZi is an orthonormal basis for L'^{Mq). If ^(xi, • ) S L^(Mo) we write 



u{xi,l) = I • dVg^. 

Jmo 

For uniqueness, let u G -f^ioco(-^) ^ solution of ()4.3p with / = 0. This 
means that for ah ip £ C^{T''^^), 

[ ui-df - 2Tdi -T^-k"^ - A^/)V' dV = 0. 
Jt 

We choose = x{xi)4>ij{x') where x e C^C^) and G C;?^(M^°*) 

with — > 0^ in H^(Mq) as j ^ cxd. Since T = R x Mq, we have 



u(xi, • dVgo {-di - 2r5i - t -k )x{xi) dxi 

R WMo 



+ / (/ OC-^x't^ijOrfV^go ) X(2;i)d2;i = 0. 

JR \JA/o / 



'Mo 

One has — Aj-'c/i^j — > Xicpi in i/~^(Mo) as j — > cx). Since n(xi, • ) G i7g(Mo) 
for a.e. xi, we have the limits as j — > 00 



u{xi, ■ )4>ij dVg,^ u{xi,l), / u{xi, ■ ){-A^f(f)ij) dVgg Xiu{xi, I), 
Mo J Mo 

which are valid for a.e. xi and for all I. Dominated convergence implies 

{-dl + 2Tdi-T'^-k'^ + Xi)u{xi,l) = inR 
for all /. By taking Fourier transforms in the xi variable we obtain 

{ei+2iT^i-T^-e + \i)u{^ul) = 0, 

with u the Fourier transform of u with respect to xi. The symbol .^^ + 2ir^i — 
— A;^ + Xi is never zero because of the condition t"^ + k"^ ^ Spec(— A^./). 

This implies that u = 0. 

Let us next show existence of solutions to (|4.3|) . We consider the case 

r > 0, the case with negative r being analogous. We start by writing (|4.4|) . 

where I is fixed, in the form 

- [{di-Tf-{Xi-k^)\u = f. (4.6) 

This equation can be factored into first order equations, where the factor- 
ization will depend on the sign of A/ — A;^. 

If Xi > k'^, then (14.60 can be written in the form 



-(Oi - r + y/Xi-k^){di - r - \/A« - k^)u = /. 



Since f £ L"^ and r / \J Xi — k'^ by assumption, Lemma 14.21 implies that 
there is a solution 

n( • , /) = -S^^^j^,S^_^^J{ -,1). 
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On the other hand, if A/ < fc^, then (j4.6p takes the form 



-(^1 - T + iy/k^-Xi){di - r - Wk^-Xi)u = f 
and since r > 1 one has a solution 

u{-,l) = -S^^^^j^S^_^^^J{ ■ , I). 

Lemma 14.21 with the trivial estimate ||f H^^s < II^IIl^, implies that 

— 6 6 



C(r+ yA7^)-i||/(-,0llL|, ^i>k', 

C{Xi-k')-^f{-,l)\\L2, Xi>k^ + AT\ 



At this point C only depends on 5. 
We write, for > 1, 

N 

un{xi,x') = '^u{xi,l)(j)i{x'). (4.7) 
1=1 

The objective is to show that as ^ cxd, u^v converges in H'^^{T) to a 
function u with ulRxaA/o = and < CT^~'^\\f\\i^2^rpy If these 

properties hold, then since u satisfies ()4.4p one has for mat 

N 

{-dl + 2Tdi -T^-k^- A,,)uN = f{xi,i)Mx'). 

1=1 

Consequently, u will be the required solution of (|4.3|) . 

Assume that t > tq > k. If0<s<2, the estimates for u show that 

oo 

Y^xt\\u{-,i)\\i, <ck'^T-' 

Z=l AKfc^ 

A:2<A,<5r2 A;>5r2 



<C7r2(-i)£||/(.,/)||2 (4.8) 
Z=l 

The last sum on the right converges since 

„ oo oo 

i=l 1=1 



Then the limit u = limiy-^oo ujy exists in LPi^{T) because of the estimate 
with s = 0. One also obtains the estimate 

For the first order derivatives, note that for fixed / 

> A;2, 



diu 



(r + ^JXl- k'^)u - 5'^_^/3^/, A, 
(r + i^Jk^ - Xi)u - S^_.^yj^J, Xi < 
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SO that • , /)||^2 < C||/( • , /)||j;^2 for T large. Also, if V^.' is the metric 

— S s 

gradient on Mq, then 

||V^'MAr||^2^(;p) = / {xiy'^^{V,j;'UN\Vx'UN)Modxi 



(Xi) '^\-Ax/UN\uN)Modxi 

R 



N 

= / y,{xiy^^Xi\uixi,l)\''dxi. 

The estimate (|4.8p with s = 1 shows that n^v converges to u in H\^{T), and 
II^IIhi^(T) — C'II/II^^Kt)- Since UN\YixdMo = 0, the same is true for u. 

Regarding the second derivatives, note that —A^iun converges in Lp^^iT) 
by ()4.8|) with s = 2, and then ||— A^^'tiH^^a ^^p) < Cr||/||j;^2(j.-). It follows 

that —/S.x'u{xi, • ) G L'^{Mq) for a.e. xi. The boundary condition for u and 
elliptic regularity imply that u{xi, • ) G H'^{Mq) for a.e. xi, and that 

l|V^/W(xi, •)llL2(Afo) < C||A^/U(xi, •)llL2(Afo)- 

Thus II V^,u||j^2 ^^j.) < Cr II / II ^2( J.). Similar estimates are true for V^'diu 
and for 

dlu = -f + {2Tdi -T^ -k'^ - A^./)m. 
This proves that ||tt||j|^2^(y) < C'r||/||2^2(2--) as required. 

It remains to consider the case of nonzero q. Denote by Gr the solution 
operator constructed above for the free case q = (we may still assume that 
c = 1), and let q be such that {xi)^^q £ L°°{T). Uu£ H^sfiC^) solves ([Ol 
with / = 0, then 

e™i(_A-A;2)e-™iu = -gu. 

Since qu G L'j(T), the uniqueness result for the free case shows that u = 
—Gr{qu) and ||it||L2 ^(j.) < Cr^-*^ ||(7u||^2(ji). If r is sufficiently large, this 

implies n = 0. For existence of a solution to (14. ip . we try u = GrV where 
V £ L'j(T) should satisfy 

{I + qGr)v = f inT. 

If r is sufficiently large then qGr is an operator on L'j(T) with norm < 1/2, 
and consequently one can take v = {I + qGr)~^ f with ||i'||j,2(2^) < 2||/||j;^2(-y). 
The required norm estimates follow from the estimates for Gr- □ 

Remark. If g G L°°{T) is compactly supported, the preceding proof shows 
that the claims in Proposition 14.11 remain true if one looks for a unique 
solution u G H^^^ q(T) instead of u G H^^^iT). 

Definition. We let Gr : Lj{T) H^^ n H\^^q{T) be the solution operator 
given in Proposition 14.11 in the case g = 0. 

In the construction of solutions, we will need to apply the operator Gr to 
functions which may not decay in xi. The proof of Proposition 14. II involves 
operators 5^ in two cases: where |Re(/Li)| is large and where Re(;u) may be 
close to 0. The latter case is problematic since good estimates may not be 
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available if there is no decay in xi. However, we only need to apply Gr to 
functions of special form: the behaviour in xi can be assumed to be like 
giAxi ^j^jgpg A > 0. The following result will be sufficient for our purposes. 

Proposition 4.3. Let A; > be fixed, let 6 > 1/2 and A > 0, and suppose 
that {xi)^^q G L°°{T). There exists tq > 1 (independent of A) such that 
whenever 

|t| > Tq and t"^ + k"^ ^ Spec(— A^/), 
then for any / = /i + /s where /i G L]{T), /s G L^^iT), and 

^xi/2( • -.x') has support in {|^i| > A} for a.e. x' G Mq, (4.9) 
there is a unique solution u G H\^q{T) of the equation 
e™i(_A -k^ + g)e~™iu = / in r. 
Further, u G H'^g{T), and the solution satisfies the estimates 



\\u\\Htg{T) < C'kl'' [II/i|Il2(T) + II/2|Il2^(T) 

with C independent of r and /i, f2- 



< s < 2, 



Proof. We follow the proof of Proposition 14. H and may assume c = 1 (this 

n+2 ri + 2 

uses that c'^ f = [c~^ fi + x/2] + /2 where x is compactly supported in 
xi, so the term in brackets is in L|(T)) and r > 0. Uniqueness is proved 
similarly as in Proposition 14. 11 and that result also gives existence if /2 = 0. 
Thus, it is enough to consider the case where f = f2- 

Assume first that q = 0. Let t > tq > k, and recall that Gr is defined by 

Grf{x^,x') = - K^,^5_,/p-^/(-,O](xi)0/(x') 

fe2<A;<5T2 

A,>5r2 

By Lemma [4. 21 the first and third sums satisfy || • ||l2(2^) < C'r~^||/||j;^2(2^) for 
any real number r. For the second sum we need to analyze the operator 
more carefully. If G R \ {0}, and if ii; G L^(R) with ii;(^) = for |^| < A, 
we have 

where mx{^) = / — fJ-)'^ and -0 is a fixed smooth function satisfying 
V' = for 1^1 < 1/2 and V = 1 for \^\ > 1. Since 

\m['\o\ < CfcA"l-^ 

we have ||5^?i'||2,2 < Cr,A 11^^11x2. 
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Using the assumption (j4.9p . the last estimate for S^, and Lemma 14.21 we 
have for k'^ < \i < 5t^ that 

~ C 

< %=^II/(-,OIIl- 

It follows that Gr maps L^iT) to L'^.iT) with norm < Cr^^. The proof that 
Gr maps into Pi HIq{T) with the right norm estimates is similar to the 
corresponding part in Proposition 14.11 

It remains to prove existence when / = /2 and q is a. potential with 
{xi)^^q G L°°{T). We seek a solution u = GrV, where v solves 

(/ + qG^)v = f. 



Here / satisfies the support condition (14. 9p . but solving this equation by 
Neumann series involves multiplication with q which breaks the support 
condition. However, we obtain a solution v = f + v if v satisfies 

{I + qGr)i = -qGrf. 

The right hand side is in L'j{T) since Gr maps / = /2 into L'^_^{T). We may 
then use the estimate in Proposition 14.11 to show that for large r there is a 
solution V e Ll{T) with ||{}||^2 < Cr-i||/||^2 . Thus, we obtain a solution 

5 —5 

to the original equation having the form 

U = Grf + GrV. 

This satisfies 111*11^^ ^(y) < GT'^~^\\f\\^2 ^^j^y □ 

Remark. As in the remark after the proof of Proposition 14. H if g G L°°[T) 
is compactly supported, the claims in Proposition 14.31 remain valid if one 
looks for a unique solution u G HI^^q{T) instead of n G H\^q{T). 



5. Norm estimates for differential forms 

The purpose in this section is to prove a counterpart of Proposition 14.31 
which applies to the Hodge Laplacian on differential forms. We will assume 
that Mq and T are as in Section HI For simplicity, we make the further 
assumptions that Mq is two dimensional (so that T has dimension 3) and 
the conformal factor satisfies c = 1. 

Let VlT = n^T e n^T e JI^T © Q^T be the graded algebra of differential 
forms. If U is in ^IT, as in Section [3] we use the vector notation 

[7 = ( i?o *5i I ^1 y (51) 

where , G O^T (j = 0, 1). 

Let —A = d5 + 5d be the Hodge Laplacian on QT, and write — A^/ = 
dx'^x' + ^x'dx' for the Hodge Laplacian on QMq. We will sometimes write 
— A^, for —Ax' acting on j-forms. Similarly to Section HI the proof of norm 
estimates for the conjugated Laplacian will require an orthonormal set of 
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eigenvectors for — A^./. In the case of 0- forms, we already used the orthonor- 
mal basis {<pi}'^i of -L^(Mo), where 

-A^>(j}i = in Mo, = on dM^. 

Here < Ai < A2 < • • • are the eigenvalues of the Laplace-Beltrami operator, 
counted with multiplicity. 

In the case of 1-forms, one needs to make a choice of boundary condi- 
tions to fix the orthonormal basis. In view of the applications to Maxwell 
equations, the relative boundary conditions (see [23l Section 5.9]) will be 
the right choice: there exists an orthonormal basis {'i/'m}m=i L^(Q^Mo) 
of real valued forms such that 

- Aj,./'0m = IJ^m^m in Mq, ttprn = ^('^x'V'm) = 0, (5.2) 

with < /ii < ;U2 < • • • the eigenvalues of — A^^' acting on 1-forms. 
We define Sobolev spaces with relative boundary values: 

HI^{^^T) = {ue H'in'^T) ; uIqt = 0} (s > 1/2), 

Hl^{n^T) = {ue H'{n^T) ■ tu = t{6u) = on dT} (s > 3/2), 

and 

i7^(QT) = {u£ H'inr) ; u has the form dM} and Rj,Sj G Hl^{n^T)}. 

We say that u is in Lj{QT) (respectively i^KOT)) if {xi^u G L'^{VlT) 
(respectively {xi)^u G H'^{^}T)). These spaces have the norms 

WuWhuut) = \\{xi)^u\\H={nT)- 

Also, u is in Hl^{nT) iff {xi)^u G Hl^i^T). 

We may now state the norm estimates and uniqueness result for graded 
forms. 

Proposition 5.1. Let k > he fixed, let 6 > 1/2 and A > 0, and suppose 
that Q : L'^_g{QT) L'j{QT) is a bounded linear operator. There exists 
To > 1 such that whenever 

|r| > To and t'^ + ^ Spec(-A°,) U Spec(-A^,), 

and whenever F = Fi + F2 where Fi G LKQT), F2 G L^_g{nT), and F2 is 
of the form F2{x) = w{xi)F2 with w a scalar function and 

supp(w;) C {|^| > A}, F2 G L°°{nT) with Vg^Fa = 0, 

then there is a unique solution U G H'^^ ^{il^T) of the equation 

e^^i(_A-A:^ + Q)e~™if/ = F in T. 

The solution satisfies the estimates 



s-l 



\^i\\r?JnT) + 11-^:: 



0<s<2, 



with C independent of r and Fi, F2. 



INVERSE PROBLEMS FOR THE ANISOTROPIC MAXWELL EQUATIONS 22 

Proof. Assume r > 0, and first consider the case where Q = 0. If is a 
0-form on T, we have already observed that 

e^i (_A - A:2)e-™iT; = {-df + 2Tdi - - - A^,)v. (5.3) 

If 7/ is a 1-form in T, we write t] = r]i dx^ + 7]' where rji = {?], dx^) and where 
rj' = r]2 dx'^ + T/3 dx"^ is a 1-form on Mq depending on the parameter xi. A 
direct computation in normal coordinates, using the identities in Section [2l 
the fact that g{xi,x') = (^ogg'^x')^^ the identity 6r] = —diiji + 5x'r]', 
implies that 

—A{7]i dx^) = {—Ar]i) dx^, 
-Arj' = -Ax'T]' - {dfrj2) dx^ - {dfrjs) dx^. 

Thus, replacing r/ by e^'^^^r] and using that Vg^^rj' = {dir]2) dx^ + (diij^) dx^, 
we obtain 

e^^i(-A - A;2)e-^^i7? = [{-df + 2Tdi - - A;^ - A^,)??i] dx^+ 

+ i-Vl^ + 2tVq, -T^-k^- A,Or/'. (5.4) 

The formulas ()5.3|) and (j5.4p give explicit expressions for the conjugated 
Helmholtz operator acting on 0-forms and 1-forms. Now * commutes with 
gT^i — fc2-jg-Txi gj]^(jg commutes with — A, so we have a corresponding 
expression for e'^^^(— A — k'^)e~'^^^ acting on graded forms written as (|5.1|) . 

Let U be as in 1^, and let F = ( \ *G° )*. Write = 

R\ dx^ + {R^y , and similarly for S^,F^,G^. The equation 

e™i(_A- A:2)e-™'iC/ = F in T 

can be written in terms of components as 

(_a2 + 2r5i -r'-e- A,,) | ^° = | , (5.5) 

(_a2 + 2Tdi -T^-k-"- A,0 I = I ^\ , (5.6) 

(_V2^ + 2rVa, - t2 - - A.,) | = | . (5.7) 

The existence and uniqueness of solutions to (15. 5p and ()5.6p follows from 
Proposition 14.31 

For the last two equations, we express {F^)' and (G^)' in terms of the 
eigenvectors ()5.2p as 

oo 

(Fi)'(xi,x') = (F^'{xum)i;^{x'), 

m=l 

oo 

(Gi)'(xi,x') = {G^'{xi,m)i^m{x'). 

m=l 
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We look for (R^)' and (S^)' in a similar form. Then (j5.7p is equivalent with 
the following equations for the partial Fourier coefficients: 

(-af+2.a, m^-^i = I (|5;(»=.''») (5.8) 

Since + A;^ ^ Spec(— A^.,), we may use the method in Propositions 14.11 
and [33] to solve 

More precisely, if {/ G H'^^^(QT) and the right hand sides in (|5.7|) are 
zero, then the relative boundary conditions imply that {—Ax'{R^y\ipm) = 
/im((-R^)'|'0m) and one obtains for all m 

Thus {{R^y\tprn) = for all m, showing that (R^Y = 0. The same argument 
applies to (S^)' ■ Existence follows by solving (jS.Sp as in Propositions 14.11 
and 14.31 and by writing (R^)' and (S^)' in terms of the Fourier coefficients. 

We have given the proof in the case Q = 0. However, the case where 
Q is bounded operator LPi^{T) — > L'j{T) is completely analogous to the 
corresponding parts of Propositions 14.11 and 14.31 □ 

6. Construction of solutions 

In this section we present a construction of complex geometrical optics 
solutions to the various Schrodinger, Dirac, and Maxwell equations which 
were introduced in Section [3l 

Let (Mo,(7o) be a simple 2-manifold, and let (Mo,go) be another simple 
2-manifold with Mq C M^^\ Write f = K x Mq, T = K x Mq, and 
rpmt _ pj^ ^ M™* for the various cylinders. Assume that T is equipped with 
the Riemannian metric 

9{xux')=(l ,° V (6.1) 



.0 go{x')^ 

The following result provides the solutions which will be used in the integral 
identity of Lemma 13.21 to recover the coefficients. Part (a) corresponds to a 
solution for the Maxwell system, and part (b) gives a solution to the Dirac 
system. We write — A^/ for the Hodge Laplacian in (Mo,go)- 

Theorem 6.1. Let {M,g) CC {T,g) be a compact manifold with boundary. 
Assume that e and /i are coefficients in M satisfying ()1.2p . (jl.Sp . (j3.ip . Let 
p be a point in Mq \ Mq, and let {r,6) be polar normal coordinates in Mq 
with center p. 

There exists tq > 1 such that for any r with 

|r| > To and t'^ + ^ Spec(-A°,) U Spec(-A;^,), 

and for any constants so,to £ the following statements hold: 

(a) For any constant A > and for any smooth function x = x(^)) there 
exists a solution to (—A — A;^ + Q)Z = in M such that one has in M 

{P -k + W)Y = 0, Y = {P + k-W^)Z, 
Y° = Y^ = 0, 
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where Z has the form 
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z 



-r(xi+jr) 



-l/4giA(a;i+ir) 



( 


So 


^ 











+ R 




to * 1 




\ 





I 





(6.2) 



and ||-R||L2(QA-f) < C\t\ ^ where C is independent of r. 

(b) There exists a solution to {P — k + W*)Y = in M of the form 



Y 



-r(xi+ir) 



-1/4 



So 

-isQ dx^ A dr 



\ 



+ R 



(6.3) 



to * 1 
\ Uq * dx^ Ndr I 

where ||-R||_L2(Q^,f) < C|t|~^ with C independent of r. 

To prove this, we begin by considering the Schrodinger equation in {T,g), 
where /c > is a constant and Q is a smooth potential with compact support 
in T™*. Following [H Section 5], we wish to construct a solution to 

(-A - A;2 + Q)Z = inT (6.4) 

by using a WKB ansatz with complex phase function, having the form 

Z = e-^P{A + R). (6.5) 

Here p = ip + is a complex weight, where (p{x) = xi is the limiting 
Carleman weight. Also, r > is a large parameter, A £ il.T is an amplitude, 
and R G H^g{Q,T) is a correction term where 6 > 1/2. 
Introduce the conjugated operators on fi'T, 

= e'Pde-^P = d-TdpA, 



5^ = e^P5e 



6 + {-ly+^T * dp A 



The conjugated Hodge Laplacian is then given by 

The next result gives explicit expressions for A,- in terms of powers of r. 
Here V is the Levi-Civita connection and Vp is the metric gradient of p. 

Lemma 6.2. If u is a 0-form or 1-form, then 



A^ii = T'^{dp, dp)u - T[2Vs7pU + (Ap)n] + A 



At-*u = * {r'^idp, dp)u - T[2V\7pU + (Ap)n] + Au} . 

Proof. The first identity for 0-forms is a straightforward computation. If u 
is a 1-form, we have 

— At-u = {d — t dpA){5u + T{dp, u)) + (5 — T * dp A *) {du — tu dp) 

= —Au + T[d{dp, u) — {5u)dp — 6{dp Au) — *dp A *du] — r'^idp, dp)u. 

The identities in Section [2] and a computation in normal coordinates show 
that 

d{dp, u) — {5u)dp — 5{dp Au) — *dp A *du = 2V\/pU + {Ap)u. 
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This proves the first identity for 1-forms. The Hodge star commutes with 
since it commutes with A, and second identity fohows. □ 

We write A = { | A^ )* where are 0-forms and 

A^,B^ are 1-forms. Using Lemma 16.21 the WKB construction for solutions 
to (j6.4p having the form (j6.5p results in the following equations in T: 



{dp,dp)=0, (6.6) 

2Vvp^^' + (Ap)AJ- = (j = 0,l), (6.7) 

2VypB^ + iAp)B^ =0 (i = 0,l), (6.8) 

e^P{-A - fc2 + Q)e-^PR = {A + k"^ - Q)A. (6.9) 

We follow the construction in [i', Section 5] and employ special coordi- 
nates to solve these equations. Considering the real and imaginary parts 
separately, the first equation (j6.6p reads 

\dtpf = \dip\^ , {dip, dip) =0. 

Recall that ip{x) = xi. Choose a point p G Mq \ Mq, and let (r, 9) be polar 
normal coordinates in (Mo,go) with center p. Then r is smooth in Mq, and 
we obtain a solution V' by setting 

ip{xi,r,6) = r. 

Note that in the {xi,r,9) coordinates one has in T 



g{xi,r,e) = I 1 
where m = \g\ is smooth. We write 



m(r, 9) 



d = -(— + i- 
2 \dxi dr 

The following result gives solutions to the transport equations ()6.7p - ()6.8p . 

Lemma 6.3. Assume the above notations. 

(1) If a is a 0-form, then 2Vvpa + (Ap)a = iff d{\g\^/'^a) = 0. 

(2) If is a 1-form, then 2Vvp?? + (Ap)?? = iff r/ = ai dx^ + ar dr + ag d9 
with 9(bl^/'ai) = 9(|<7|i/V) = d{\g\-"^ae) = 0. 

Proof. We have 

p = xi+ir, Vp = 2d, /S.p = d \og\g\. 

The equation for 0-forms is 49a -|- {d log l^l) a = 0, which proves the first 
part. For the second part, the form of g shows that 

pi ^ Q r' = I I^D' k = l = 9, 

-^^ ' rfe I otherwise. 
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Consequently Vg^ dx^ = for all j and Vg,, dx^ = Vq^ dr = 0, Vg^ dO = 
— ^9r.(log \g\)d9. The result follows by noting that 

VvpCo^i dx^ + ttrdr + ag dO) 

= {2dai) dx^ + {2dar) dr + {2dag - ^^^(log \g\)ae) dO 

where |9r(log \g\) = 9(log \g\). □ 

We are now ready to give the construction of complex geometrical optics 
solutions to the Schrodinger equation. 

Proposition 6.4. Let (Mq, 50) CC (Mo,(7o) be two simple 2-manifolds, and 
consider the cylinders T = R x Mq and T = R x Mq equipped with the 
metric g given by ()6.ip . Let > be a constant, let b > 1/2, and let Q be 
a bounded linear operator L'^^{QT) L'j{^T). There exists tq > 1 with 
the following property: if 

|r| > To and t'^ + ^ Spec(-A°/) U Spec(-A;^,), 

and if p, A, ao, oi, 60, 61 are any parameters such that 

p is a point in Mq \ Mq and A > is a constant, 
(r, 6) are polar normal coordinates in Mq with center p, 
a;, 6; are smooth functions in T of the form e'^^^'"w{x'), 
{di + idr)ai = {di + idr)bi = in T, 

then the equation (—A — A;^ + Q)Z = in T has a unique solution 

(6.10) 

where R G H'^^^{QT). The remainder R satisfies ||i?||j;^2 (q^-^ < C|r|~^ 
with C independent of r. 

Proof. Take p = xi + ir and = \g\~^^^aQ, = \g\~^^^aidx^ , and also 
= \g\~^^'^bo, = \g\~^^^bi dx^. It follows from the discussion above that 
equations (j6.6|) - (j6.8|) are satisfied, and that Z solves (—A — k'^ + Q)Z = 
iff 

e'^f{-A-k^ + Q)e-^PR = F (6.11) 

where F = (A + A:^ - Q)A. 

The form of ai and hi implies that A £ L'i^{QT), so QA £ Lj{^T), 
and also that (A - k^)A = e'^''^F2{x') G L'i^i^T) where F2 G L°°{nT), 
'^8x^2 = 0. The latter fact follows from the formula for A acting on VtT 
computed in the proof of Proposition 15.11 Thus we have a decomposition 
F = Fi -\- F2 as in Proposition 15.11 and that result shows that there is a 
unique solution R £ H'^^j^{ilT) to (j6.1ip with the required estimate if |r| 
is large and outside a discrete set. □ 





1 ao \ 






hi * dx^ 


+ R 


60 * 1 
\ ai dx^ ) 
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We now prove the main result on complex geometrical optics solutions. 
For part (a) we need to use the uniqueness of the solutions above to conclude 
that = = 0. Part (b) is in fact much easier since it is enough to 
construct solutions to a Dirac equation in M without worrying about the 
vanishing of scalar parts. 

Proof of Theorem. \6. l]f a) . Assume the conditions in Theorem 16.11 and ex- 
tend e and fi smoothly as constants into T. Then Q satisfies the assumption 
in Proposition 16.41 and that result guarantees the existence of a solution Z 
to (-A -k^ + Q)Z = in T of the form (IHIH with R G H'i^^^{T) and 
^{nx) < C\t\-^. Setting Y = {P + k- W*)Z, Lemma O shows that 
Y solves {P-k + W)Y = in T. 

The main point is to show that Y^ = = 0. For this we use an idea 
appearing in [19j . By Lemma [3. II we have (—A — k'^-'r Q')Y = 0. Looking at 
the 0-form and 3- form parts and using the special form of Q' , the equation 
decouples and we obtain the following equations in T: 

(-A - A;2 + gO)yO = 0, 
(-A -k^ + q^) * y3 = 0, 

where and are smooth potentials with compact support in M™*, 

q^ = k^ — — 2^0? + -(dcK, da). 
Now, writing p = xi + ir and Z = e~'^^{A + R), Y^ has the form 

yo = ((P + A; - W')Zf = e'^P {{^-\^^ + k-W'^{A + R)^ 
= e-^P {yo + ro) . 
Here we have written 

yO = -^SrA' +kA'>, 
I 

ro = -(W\A + R)f - -6rR^ + kR^. 

Since R G H^^^^QT) we see that ro G H\^{T) and rolar = 0. 
We will choose AP , A^ so that 

yo = 0. (6.12) 
Then e~*'^''ro will be a solution in H^^^{T) of the equation 

e-^i (_ A -k^ + gO)e-™i (e-*"Vo) = in T, 

and the uniqueness part in Proposition 14.11 will show that ro = 0, so also 
Y^ = 0, if |t| is sufficiently large. To obtain (|6.12p we make the choices 
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= \g\~^/^ao, = \g\~^/^ai dx^ where 
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Ol 



ao 



ikso 



xi0), 



—5r{ai dx^] 
ik 



(6.13) 
(6.14) 



This is consistent with Proposition 16.41 since {di + idr)ai = and since 

1 T — iX 

ao = —{-diai + rai) = — - — oi (6.15) 



so also [di + idr)ao = 0, and both oq and oi are of the form e'^^^^w{x'). Now 
(|6.12|) holds because 

yo = --6rA^ + kA^ = --6r{\g\-^/^ai dx^) + k\g\-^'^^5r{ai dx^) = 0. 
i i ik 

We have established that = given the choices (I6.13p ~ (l6.14|) . A 
similar computation for Y^, with 



bi 
bo 



ikto 



xio), 



ik 



5r{bi dx^] 



shows that Y^ = 0. Finally, we note that Z is of the form 

Z = e-^^il^r^/^ ( ao *bi dx^ \ *bo ai dx^ f + R] 

where oq = soe*^(^i+*^)x(f ) + 0(|r|-^) by ^JEh . and oi = 0{\t\~^). We 
have written 0(|r|~^) for quantities whose LF'^M) norm is < C|t|^^. Similar 
expressions are true for 6o and bi. This shows that Z has the required form 
^ with \\R\\L2(nM) < C\t\'\ 

(Note that aj and bj are mildly r-dependent, but their W^'°°{T) norms 
are bounded uniformly in r which implies that final constant C does not 
depend on r). □ 

Proof of Theorem \6.1]( h). Again, assume the conditions in Theorem 16.11 and 
extend e and /U smoothly as constants into T. Let Q be the potential in 
Lemma [3. 1[ We look for a solution to (—A — k'^ + Q)Z = in M of the form 



Z 



3-r(xi+jr) 



-1/4 





to * dx^ 



\ So dx^ ) 



+ R 



Write p = xi + ir and A = \g\ ( to * d'^^ \ sodx^ )*. Following 
the WKB construction, it is enough to solve 

e^P(-A - A:^ + Q){e-^PR) = {A + k^ - Q)A in M. 

Define F G L^(Or) with F = e-*^'^(A + A:^ - Q)A in M and F = in 
T \ M, and let e~'^'^^R be a solution provided by Proposition 15.11 of the 
equation e™'^ (— A — /c^+Q)(e~'^^^ [e~^'^^R]) = F in T. This gives the required 
solution Z in M satisfying ||-R||L2(f2jvf) ^ C'|r|~^ and ||-R||_f/i(nAf) ^ C. 
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We set 



Y = -{P + k-W)Z. 



By Lemma O this satisfies {P - k + W*)Y = in M, and 
Y = e-^P (-^dr --^5r + k-w\ {t'^A + t^^R) 



-TP 



\a\ 



-1/4 



-1/4 



dp A (so dx^ + to * dx^) 
*dpA *{so dx^ - to* dx^) + 0(|t|"^) 



-TP 



( 



-1/4 



ISO 



So dx^ A dr 



ito * 1 
\^ -to * dx^ l\dr j 



+ 0(|t|-i1 



Here 0(|t| ^) denotes a quantity whose l?(yiM) norm is < C\t\ ^. The 
result follows upon replacing so by —iso and to by — ito- □ 



7. Recovering the coefficients 

We shall use the complex geometrical optics solutions constructed in The- 
orem [6TT] to prove Theorem ll.il The first step is a reduction to the case where 
the conformal factor in the metric is equal to one. We write A = ^g^e.p. for 
the admittance map in {M,g) with coefficients e and p. 

Lemma 7.1. Let {M,g) be a compact Riemannian 3-manifold with smooth 
boundary, and let e and fi satisfy (|1.2p - ()1.4p . If c is any smooth positive 
function on M, then 



A 



Proof. Follows by noting that *cgU = c^/^^^ *g u for a /c-form n, so that 
a pair [E, H) satisfies (jl.ip with metric eg and coefficients e and /x iff it 
satisfies (jl.ip with metric g and coefficients c^^'^e and c^/'^p,. □ 

Proof of Theorem li.il According to Lemma 13.21 we may assume that p.7p 
holds and the identity ()3.8|) is valid. By the definition of admissible mani- 
folds, there are global coordinates x = such that g has the form 



g{x) = c{x) 



1 

go{x'] 



If A„ 



A„ 



J, then also A i -/ 



^g,ei,p.i — ^^g,e2,p.2-. ^^^^^^ '^c-^g,c^/'^ei,c^/'^p.i " ^c~^g,c^/'^e2,c^/'^p.2 

Lemma l7. II This shows that we may also assume c = 1. 

By Theorem 16.11 if r is outside a discrete set and |t| is sufficiently large, 
and if x(^) is a smooth function and A > and sq and to are real numbers, 
then there exist Zi and Y2 in H'^(QM) satisfying the conditions in Lemma 
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and having the form 
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-T(x\+ir) 



^T{x\—ir) 



( 



-l/4giA(xi+ir) 



So 




to * 1 




\9\ 



-1/4 



So 

iso dx^ A dr 



\ 



+ R2 



to * 1 
\ —ito * dx^ f\dr J 

where ||L2(r2M) ^ C'|t|^^ with C independent of r. In the second solu- 
tion, we used — r instead of r as the solution of the eikonal equation. 
By Lemma 13.21 these solutions satisfy the identity 



{{Qi-Q2)Zi,Y2)dV = {). 



M 



Letting r ^ cxd outside the discrete set and using the estimates for we 
obtain in terms of the x = (xi, r, 9) coordinates that 




Q2 



/ So \ / So \ 

—isQ dx^ A dr 



to * 1 io * 1 

\ / \ ito * dx'^ A dr J 

Let Qa and qp be the elements of Q1 — Q2, interpreted as a 8 x 8 matrix, which 
correspond to the (1, l)th and (5, 5)th elements, respectively. By Lemma [3. II 



^A(qi - 02) + ^{dai,dai) - ^ 
1 . _ „ . 1 , , 1 



9/3 



{da2,da2) -u;^(ei/xi -£2/^2)! 



2-A(/3i - /?2) + ^{dPudPi) - ^{dP2,dP2) 

With the two choices (so,to) = (1,0) and (so,to) = (0,1), the special form 
of Qi and Q2 in Lemma l3. II shows that we obtain the two identities 

2*^(^1+*'') X(^)9«(x)rfx 



M ' 



0, 
0. 



We extend Qa and Qjs to be zero in T \ M, where T DD M is as in the 
definition of admissible manifolds. Then the integrals above may be taken 
over T = R X Mq. Varying x(0), it follows that for all 6 we have 



-Xr 



z'^''^qaixi,r,e)dxi 



dr = 



and similarly for Q/^. Now, since {r,6) are polar normal coordinates in Mq, 
the curves r (r, 6) are geodesies in Mq. Denoting the expression in brack- 
ets by fa{r,9) and varying the point p in Theorem 16.11 and varying 9, we 
obtain that 



/a(7(?'))exp 



A ds 



dr = 



for all geodesies 7 in Mq which begin and end at points of OMq. This shows 
the vanishing of the geodesic ray transform of the function fa with constant 
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attenuation —A. For more details we refer to [U Section 7]. In particular, 
the injectivity result given by Theorem 7.1 in [1] implies that = for all 
positive A which are sufficiently small. Thus 



/oo 
-oo 



for such A and for all r and 6. Since Qa is compactly supported in xi, the 
Paley- Wiener theorem shows that = in M. We obtain = in M by 
the exact same argument. 

We have arrived at the following two equations in M: 

-^A{ai - 02) - l{d{ai + 02), d{ai - 02)) + '^^(ei^i - £2/^2) = 0, 
-iA(/3i - 132) - j{d{(3i + (32),d{pi - P2)) + - ^2/^2) = 0. 

Let u = (£1/62)^''^ and v = {^1/ ^2)^^"^ ■ Then ^(ai — 02) = log u, and the 
equations become 

-A(log u) - {eie2)-^''^{d{eie2Yl'^,d{\og u)) + Lo'^{ei^ii - E2ti2) = 0, 
-A(log v) - (^l^2)"^/^(^i(/UlAi2)^/^d(log v)) + Lo'^{eiiii - e2ti2) = 0. 

Multiplying the first equation by (£162)"^''^ and the second by {^11^2)^^^, and 
using that 6{aVw) = —aAw — {da,dw), we see that u and v satisfy the 
semilinear elliptic system 

6{e2du) + Lu'^elfi2iu'^v'^ - '^)u = 0, 
5{iX2dv) + uJ'^e2^l2{u^v'^ — l)v = 0. 

The condition (|3.7p ensures that one has n = 1 and v = 1 near dM . Also, 
the above equations imply that the pair {u,v) = (1, 1) is a solution of the 
semilinear system in all of M. By Theorem IB . 1 1 uniq ue continuation holds 
for this system, and we obtain n = 1 and u = 1 in M. This proves that 
El = £2 and fii = ^2 in M as required. □ 

We now prove Theorem 11.21 The treatment below follows [12j. Let 
Q, C r3 be a bounded open set with smooth boundary, and let e and 
be symmetric positive definite (1, l)-tensors on ^1. We equip with the 
Euclidean metric e. The Maxwell equations (jl.Sp can be written as 



curie (£;) = iuJ^H, 
curie (if) = —iueE. 
Here 

curle(X) = {*edX^f 
with the flat and sharp operators taken with respect to e. 

For the vector fields E = {Ei,E2,E3) and H = {Hi,H2,H3), let E = 
& = Ej dx^ and H = = Hj dx^ be the corresponding 1-forms. To write 
(|7.1|) in a form similar to p.l|) . it is enough to find Riemannian metrics 
and so that 



* 



These conditions will be satisfied if we choose in local coordinates 
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Then (j7.ip is equivalent with 

*g^dH = —iujE. ^ ' 

We now use the assumption that e and [i are in the same conformal class, 
so that ^ = o?e for some smooth positive function a on fi. This allows to 
define a metric (7 on 17 by 

g = a^ge = a~'^g^. 
Since *cgU = c^^/^ *g u for a 2-form u, (j7.3|) is equivalent with 



Kgdi? = iivaH. 
'9 



This is of the form (II. ip . and further the tangential boundary condition 
tE = f is of the same form as ()1.6p . Since ()7.4|) is equivalent with ()1.5p . 
Theorem lA.ll implies that for uj outside a discrete set of frequencies the 
system (|1.5p ~ (|1.6p is uniquely solvable for a given boundary value /. The 
admittance map A for (jl.Sp reduces to the map 

Given this reduction, it is easy to prove the second main result of the 
paper. 

Proof of Theorem \1.2l Upon interpreting e and fi as (0, 2)-tensors by raising 
one index with respect to the Euclidean metric, the condition (|7.2p implies 
that 

e = det{e)g~^, fi = det{fi)g~^ . 

From the assumption in the theorem, we know that there is an admissible 
metric g and smooth positive functions cj ,Cj on for which 

9e, = c]g, gf,^ = c]g. 

Using these formulas in (17. 3p . the Maxwell equations in Vt with coefficients 
Ej and are equivalent with 

*gdE = iujCjH, 
*gdH = —luJCjE. 

Since the admittance maps for the Maxwell equations in 17 coincide, it follows 
that Ag^ci,ci = ^g,c2,c2- From Theorem 11.11 we obtain ci = C2 and ci = C2, 
which implies that gi.-^ = g^^ and g^^ = g^^ ■ By (I7.2p 

1 _ 1 1 _ 1 

det(ei) det(e2) ^' det(^i)^^ det(/i2)^^ 

Taking determinants gives that det(ei) = det(e2) and det(/ii) = det(/i2)- 
Consequently ei = £2 and = ^2- D 

Remark. Note that in the setting of (jl.ip . if e and fi are real valued, then 
a conformal scaling of the metric would reduce (jl.ip to a system of the form 
(j7.4p . Therefore, in Sections [3] and [6] it would be enough to consider the case 
where fi = e~^, which would simplify some of the arguments slightly. 
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Appendix A. Wellposedness theory 

Let {M,g) be a compact oriented Riemannian 3-manifold with smooth 
boundary dM. Consider the Maxwell equations 

*dE = iujfiH in M, 



*dH = -iujeE in M, ^^'^^ 

with the tangential boundary condition 

tE = f on dM. (A.2) 

Here we assume that e and ^ are complex functions in C^{M) whose 
real parts are positive in M, and w is a complex number. To describe the 
boundary condition in more detail, we introduce the Div-spaces 

Hf)i^{M) = {u£ H'Q}{M) ; Div(t'u) G H'-^/'^{dM)}, 

THl,-„{dM) = {/ G H'n\dM) ; Div(/) e H'idM)}. 

There are Hilbert spaces with norms 

ll^ll//£,i^(M) = \\u\\H^n^{Ai) + l|Div(tu)||^s-i/2(g^), 

\\f\\THf^^^{dM) = WfllH^Q^dM) + l|Div(/)||^s(5M). 

It is easy to see that t{Hl^^^{M)) = TH^^^''^{dM) for s > 1/2. 

Theorem A.l. Let e,/i G C^{M), /c > 2, be functions with positive real 
parts. There is a discrete subset S of C such that if lv is outside this set, 
then one has a unique solution {E, H) G H^-„{M) x H^-^^{M) of ^0) ^ (11:2]) 

given any / G TH^.J^^idM). The solution satisfies 



with C independent of /. In particular, if e, G C°°{M) and / G ^^{dM), 
then there is a unique solution {E,H) G r2^(M) x 0^(M). 

The existence of a solution will be proved by the well-known variational 
method as in [3], [13]. We proceed to describe this method. The first step is 
to solve for H in the first line of (lA.ip and to substitute this on the second 
line, which leads to the second order equation 

6{fi-^dE) - u'^eE = 0. 

However, this equation does not imply the divergence condition 5{£E) = 
which is necessary for solutions of (jA.ip . To make up for this, we consider 
the modified equation 

6{fi-^dE) + se-^d5{eE) - u'^eE = 

where s is a positive real number. The condition 6{£E) = will follow later 
from the equation 

s6{e'^d[6{eE)]) - uj^5{eE) = 0, 

which is obtained by applying 6 to the earlier equation. 

To connect the present situation to boundary value problems for the 
Hodge Laplacian, we write e = eE and note that e should satisfy 

6{fi-^d{e~^e)) + se-^Se - uj'^e = 0. (A.3) 



INVERSE PROBLEMS FOR THE ANISOTROPIC MAXWELL EQUATIONS 34 

Taking the LP' inner product of this with ee for a 1-form e, and assuming 
the relative boundary conditions (see [231 Section 5.9]) 

te = te = and t{6e) = t{6e) = 

or the absolute boundary conditions 

t{*e) = t{*e) = and t{6 * e) = t{5 * e) = 0, 

we end up with the following bilinear form for solving the Maxwell equations. 

Definition. If e, e G Hln^{M) (where b = Ror b = A), we define 

B{e, e) = {fi-^d{e-^e)\d{ee)) + s{6e\6e). 

Here we have used the spaces 

H\9}{M) = {n G H^n^{M) ; tu = 0}, 

H\^^{M) = {n G H^n^{M) ; t{*u) = 0}. 

If A; > 2 we define H'lfl}{M) = {u G H''n^{M); tu = t{5u) = 0} and 
H\il^{M) = {n G H^n^{M) ; t{*u) = t{6 * u) = 0}. 

One defines weak solutions of ()A.3P in the usual way. The main point is 
the following solvability result. 

Proposition A. 2. Let e and /j, be functions in C^(M) with positive real 
parts, and let s be a positive real number. There is a discrete set in C 
such that if lu is outside this set, then for any / G {H^il^{M)y the equation 

6{fi-^d{e-^e)) + se'^dde -uj^e = f (A.4) 

has a unique solution e G HI^Q^{M) {b = R or A). One has 

\\e\\m<C\\f\\^Hiy. 

Proof. Clearly i? is a sesquilinear form on H^i^^{M) with 

\B{e,e)\ < C\\e\\fji\\e\\fji. 

We may write B{e, e) = i?o(e, e) + -Bi(e, e) where 

Bo{e,e) = {fi~^de\de) + s{5e\6e), 

and i?i is a sesquilinear form such that |Si(e,e)| < C||e||2,2 ||e||j|/i . It follows 
that 

ReB{e,e) > c||de||^2 + s||5e||^2 — C||e||2,2 ||e||j|;-i . 

We now invoke a Poincare inequality for 1-forms with relative or absolute 
boundary values: by [23l Section 5.9] one has 

ll^ll^i < C{\\u\\l, + \\du\\l, + \\6u\\l,), u G Hln\M). 

It follows that 

Rei?(e, e) > c||e||^i — C||e||^2. 

We have proved that for some Cq > 0, the sesquilinear form B + d is 
bounded and coercive on H^Q^{M) for any d G L°°{M) with Re(d) > 
Cq. By the Lax-Milgram theorem there is a bounded linear operator T : 
{H^^l^y — > Hi,^^ which maps / to the unique solution e of (|A.4p where 
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— is replaced by the constant Ci = Co/mina;gM Re(e). Now, e solves 
(II3D iff 

{I-{Ci+LO^)T)e = Tf. 

The last equation has a unique solution iff either (Ci + ^ Spec(T) or 

cj^ = — Ci. The operator T : H^Q,^ ^1^^ compact by the compact 
embedding HIQ} — > L^il^, and ^ Spec(r), so Spec(T) is discrete. Then 
the set 

= {(J G C \ {±i^/Ci} ; {Ci + e Spec(T)} 

is also discrete and ()A.4p is uniquely solvable for any w ^ C \ Eg. □ 

Given the last result, higher order regularity for solutions follows in a 
similar way as for the Hodge Laplacian (for more details see \1?>\ Proposition 
9.7] and the results mentioned there). 

Proposition A. 3. Let e and /i be functions in C^{M), k >2, with positive 
real parts, and let s > 0. If a; ^ S^, then for any / G H''~'^fl^{M) the 
equation (IA.4|) has a unique solution e G H^Q^{M) {b = R oi A), and 

\\e\\H,<C\\f\\H.-2. 

Further, if u; is any complex number and if e G H^Q^(M) solves (|A.4p for 
some / G H^-'^n^{M), then e G H^9}{M). 

Finally, we connect the above discussion to the Maxwell system and prove 
the main result. 

Proof of Theorem \A.l[ We take S to be the set Si in Proposition IA.21 and 
assume that u; ^ S. As a technical preparation, we choose s > so that 
Lo'^/s is not an eigenvalue of the operator u ^ 6{e~^du) defined on Hq{M). 
We then have u; ^ S<j, which may be seen as follows: if e G Hj^Q^{M) 
satisfies (|A.3p . then e G Hj^Q^{M) by Proposition IA.31 and applying 5 to 
both sides of (\A.3h shows that u = 6e is a solution in Hq(M) of 

s6{e^^du) - uP'u = 0. 

By the choice of s we have u = 0, which implies that e satisfies ()A.3p with 
s = 1. Then e = by the assumption uj ^ Si, and therefore oj ^ S^. 

For uniqueness, if {E,H) G H}^-„{M) x H}^-„{M) solve (fXT]l - (|X2l) with 
/ = 0, then one has 

5{iJL-^dE) = Lo'^eE, 
6{eE) = 0. 

It follows that e = e^; is in H}fL^{M) and 

6{^i'^d{e~^e)) + e'^dSe - u?e = 0. 

Proposition IA.2I shows that e = 0, which implies E = H = 0. 

Let us proceed to prove existence of solutions. Given / G TH^^^^'^{dM), 
choose G H^^^{M) with tE^ = f and t5{eE^) = 0. A computation in 
boundary normal coordinates shows that the extension map f ^ E^ can be 
taken to be bounded and linear TH^^''^{dM) H^^^{M). 
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We let e G H'^Q^[M) be the solution, given by Proposition IA.31 of 

6{n-^d{e-^e)) + se~^d6e -Lo^e = F 

with F = -5{^x-^dE^) - se~^d5{eE^) +lo'^£E° £ H''-^n^{M). Now define 
E = e~^e + E^ and H = J— * dE. With these conventions, E satisfies the 
equation 

6{ii-^dE) + se-^d6{eE) - uj^eE = 0. (A.5) 
We now claim that 

6{eE) = 0. (A.6) 
In fact, by taking 6 of both sides of (|A.5p . the function u = 5{eE) G Hq{M) 
satisfies s6{e^^du) — uj'^u = 0, showing that u = by the choice of s. 

The first equation in (jA.ip is satisfied by definition, and also the second 
equation is valid since by ([Oil and (IXBll 

*dH = —6(n'^dE) = -—e'^d6(eE) - iueE = -itoeE. 

ILO lUJ 

The 1-form E is in H^^^{M) and tE = f. The 1-form H is initially in 
H^^^il^{M) by definition. However, a similar argument which was used to 
establish (IA.3|) shows that h = fiH satisfies the second order equation 

6{e-^d{fi-^h)) + fi-^d6h - uj^h = 0. 

Also, a computation in boundary normal coordinates gives the following 
boundary conditions for h: 

tUh) = —tidE) = — Div(/) dS, 

ILO ILO 

t{5 *h) = t{*d{nH)) = t{*dfi AH)- iujfxef. 

Since / S TH^^^^'^ and H £ H^^^^ one can check by a computation in 
boundary normal coordinates that there exists hP G H^Q}{M) for which 
h = h — hP is in H^^Q^{M). Now h satisfies the equation 

5{£-^d{^x-^h)) + ^i-^d5h - uj'^h = F 

for some F G H^~'^Q}{M). Elliptic regularity (as in Proposition I A . 3p implies 
that h G H^n^{M) which is then true for H too. We have H E H^^^{M) 
because 

Biv{tH) = -{u,*dH)\QM = iu;e{iy,E)\9M G H^~^/\dM). 

□ 

Remarks. 1. If e, /i G C^(M), then the conclusion of Theorem lA.il is 
valid also for k = 1. This follows from the above argument upon 

1/2 

approximating / G TH-^-^^ by smooth tangential fields. 
2. Theorem I A . 1 1 consider s the case where e and /x are independent of w. 
In applications the coefficients are often w-dependent, for instance in 
lossy materials one writes e = Re(e) + ia/u; with Re(e) > 0, cr > 0, 
and > independent of w. In the last case, wellposedness in 
Euclidean domains was shown in ^Oj. Of course, Theorems 11.11 and 
11.21 are valid whenever the admittance map is well defined (this is 
the content of assumption (II. 4p ). 
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Appendix B. Unique continuation 

This section contains a unique continuation result for principally diagonal 
systems required in the final recovery of coefficients. The result is well known 
and follows from standard scalar Carleman estimates, but since we could not 
find a proper reference a proof is included here. 

Theorem B.l. Let {M,g) be a compact connected Riemannian manifold 
with boundary, and let ar,fir be Lipschitz continuous functions in M with 
positive real parts (r = 1, . . . , N). Consider the operators 

PrU = — S{l3rdu), 

Or 

P = diag(Pi,...,P7v). 
Let r be an open subset of dM. If u G H^{M)^ satisfies 

\Pu{x)\ < C{\u{x)\ + |Vn(x)|) for a.e. xe M, 
u\r = d^u\r = 0, 

then n = in M. 

More generally, strong unique continuation holds in this setting. We will 
deduce Theorem IB. II from the next result which is stated in R". 

Theorem B.2. Let B he a ball in R" with center xq, let {g^^)"^ i^^i be a 
Lipschitz continuous symmetric positive definite matrix in and let ar,Pr 
be Lipschitz continuous functions in B with positive real parts. Consider 
the operators 

a.r 

P = diag{Pu ■■■,Pn). 
liuG H^{B)^ satisfies for aU > 

\Pu{x)\ < C(\u(x)\ + |Vn(a;)|) for a.e. x e B, 

lim / \u{x)\'^ dx = 0, 

'•^O JB(xQ,r) 

then u = m B. 

Proof of Theorem \B.ll If u is as in Theorem IB.lt we fix a point on T and 
take M to be a manifold obtained by enlarging M slightly near this point. 
Extending u by zero to M and extending and Pr as Lipschitz functions, 
we see that |Pn| < C{\u\ + |Vu|) a.e. in M and n = in some open subset. 
Working in local coordinates and using Theorem IB. 21 with a connectedness 
argument proves the result. □ 

To prove Theorem IB.2( note that we can assume ar = 1, and by dif- 
ferentiation that f3r = I- Letting Lu = dxjig-'^'dxf.u), this implies that 
Pi = . . . = Pn = L. We may also assume that xq = and B = B{0, 1), 
and that u is real valued. The result is a consequence of the following scalar 
Carleman estimate. 
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Proposition B.3. Let A > be such that A^^l^p < g^^^j^k < A|CP in B 
and X^"fc=il5'"''^(^) ~ 9^^{y)\ ^ -^1^ ~ ul x,y ^ B. There exists < 6 < 1 
and M > 0, only depending on n and A, and a function w satisfying 

\x\/M < w{x) < M\x\ in B, 

such that for all a> M and all u G C^{B{0,6) \ {0}) we have 



This estimate is proved in [51 Theorem 2.1], and is also contained in 
[S] with slightly different hypotheses. Theorem 2.1 in [5J is given in the 
parabolic setting, but the result above follows by applying the estimate 
in [3] to v{x,t) = 6{t)u{x) where is a cutoff function. One then notes 
that when the coefficients are independent of t, the w{x, t) constructed in 
[5] depends only on rr, so the inequality (2.1) in [5] yields the lemma by 
integrating in t and absorbing the extra term on the right by making a even 
larger. 

As an immediate corollary to Proposition IB . 31 for u G C^(-B(0,(5)\{0})^ 
we have that 



Once we have the last estimate, after the initial reductions. Theorem IB. 21 
follows immediately using the standard Carleman method. 
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